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Abstract 

In the prospect of diffractive Higgs production at the LHC collider, we give an extensive study 
of Higgs boson, dijet, diphoton and dilepton production at hadronic colliders via diffraction at both 
hadron vertices. Our model, based on non factorizable Pomeron exchange, describes well the observed 
dijet rate observed at Tevatron Run I. Taking the absolute normalization from data, our predictions 
are given for diffractive processes at Tevatron and LHC. Stringent tests of our model and of its 
parameters using data being taken now at Tevatron Run II are suggested. These measurements will 
also allow to discriminate between various models and finally to give precise predictions on diffractive 
Higgs boson production cross-section at the LHC. 

1 Double Diffractive Hard Production through non-factorizable 
Vacuum Exchange 

The discovery of the Higgs boson is one of the main goals of searches at the present and next hadronic 
colliders, the Tevatron and the LHC. The standard, non-diffractive production mechanisms are being 
studied extensively. The main decay modes of a low-mass Higgs boson are 6-quark pairs and r-leptons, 
which are difficult to extract from the standard model background processes. A promising channel is the 
77 decay mode; however, due to the small branching fraction (about 10~ 3 ), a luminosity of the order of 
50 fb _1 is needed to establish the existence of the Higgs boson in the mass range below 135 GeV. It is 
thus important to investigate other, complementary ways to produce and detect the Higgs boson in this 
mass range. 

One promising production mode, the exclusive double diffractive production, was proposed some time 
ago in Ref. In this case, the Higgs boson is diffractively produced in the central region resulting in 

a final state composed of the two protons scattered at very small angles and detected in the roman pot 
detectors, the decay products of the Higgs boson in the main detector, and nothing else. It is thus a very 
clean signal. The kinematic constraints coming from the proton detection in the roman pot detectors 
allow a very precise determination of the Higgs boson mass [3] , hence improving the signal to background 
ratio. Contrary to non-diffractive production, the main Higgs boson decay modes, like bb or rr are 
thus promising channels. However, the exclusive cross-sections may be very low, and thus put strong 
limitations to the potentialities of double diffractive production. 

*CERN, CH-1211, Geneva 23, Switzerland and CEA/DSM/DAPNIA/SPP, CE-Saclay, F-91191 Gif-sur-Yvette Cedex, 
France 

tCEA/DSM/SPhT, Unite de recherche associee au CNRS, CE-Saclay, F-91191 Gif-sur-Yvette Cedex, France 
tCEA/DSM/DAPNIA/SPP, F-91191 Gif-sur-Yvette Cedex, France 



1 



Recently, we have studied the possibility of producing the Higgs boson together with other particles, 
i.e. inclusive diffractive production 0]. The expected cross-sections are increased compared to exclusive 
production, and the model can be "calibrated" using the diffractive dijet production measured by the CDF 
Collaboration at Tevatron run I [5]. The experimental result, and in particular the dijet mass fraction 
spectrum, shows that hadronic activity in the central region (coming from hard QCD radiation, and from 
soft Pomeron remnants) needs to be accounted for, in order to describe these data. We have proposed a 
model for double diffractive production of heavy "objects" , based on a non-factorizable Pomeron model 
and able to describe the observed features of dijet production data in a qualitative way. Normalizing 
our raw predictions to the CDF measurement allows to make quantitative predictions for the Tevatron 
and the LHC, given specific assumptions for the model parameters. We have also shown [S] that it is 
possible to reconstruct precisely the mass of the Higgs boson if both protons in the final state can be 
detected with roman pot detectors and if the Pomeron remnants can be measured in the forward region 
with sufficient resolution. 

At the Tevatron collider, the Higgs boson production cross-section is severely limited by the small 
available phase space. At the LHC, since the beam energy is much higher, diffractive Higgs boson 
production has larger cross-section, and might thus be an interesting channel, as we shall discuss further 
on. On the other hand, even if double diffractive Higgs boson production at the Tevatron is probably too 
small by itself, it is however possible to verify and constrain models of double diffraction at the Tevatron 
Run II, through the study of difcrmion production. 

In this paper, we wish to provide an extensive study of our model, including predictions for dijet, 
diphoton and dilepton production and their ratios for both the Tevatron run II and the LHC. We will 
discuss possible values of its characteristic parameters, and propose means to discriminate between the 
existing models at the Tevatron. In turn, this study will help to formulate more precise predictions for 
double diffractive Higgs boson production at the LHC in the near future. 

The plan of our study is the following. In Section |3 we will recall and discuss the original exclusive 
model of Bialas and Landshoff PP, and our extension to inclusive production. In Section [3] predictions 
are given for dijet, diphoton, dilepton and Higgs boson cross-sections. Dijet results for the Tevatron run 
I, using a fast simulation of detector effects, are in good agreement with available data and are used to 
normalize our predictions. A discussion of Pomeron remnants and their possible detection is also given. 
In Section^] ways to verify our model and determine its parameters using Tevatron data are proposed. 
The possibility of distinguishing our model from the existing models of double diffraction is discussed in 
Section [SJ Section [H] summarizes the interest of diffractive Higgs boson production compared to standard 
production. Finally, Section [3 gives conclusions and an outlook on the promising future studies on hard 
diffraction production at the Tevatron and the LHC. 

2 Inclusive vs. exclusive production 
2.1 Exclusive production 

Let us first introduce the original model of ^] describing exclusive Higgs boson and qq production in double 
diffractive production (noted DPE 1 in the following). This process is depicted in Fig. ^ "Exclusive" 
case. 

keep the now standard notation DPE (i.e. Double Pomeron Exchange) but the model we describe, together with its 
extension to inclusive production, considers a non-factorizable soft pomeron exchange with one common exchanged gluon, 
while a factorizable double Pomeron mechanism implies two different pairs of gluons (see Section HI. 
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In PP, the diffractive mechanism is based on two-gluon exchange between the two incoming protons. 
The soft Pomeron is seen as a pair of gluons non-perturbatively coupled to the proton. One of the gluons 
is then coupled perturbatively to the hard process (either the Higgs boson, or the qq pair, see Fig. 
while the other one plays the role of a soft screening of color, allowing for diffraction to occur. This soft 
character requires the phenomenological introduction of a distinctive non-perturbative gluon propagator 
[2] whose parameters are constrained by the description of total cross-sections using the same formalism. 

The hard gluons, carrying all remaining momentum (xf = x g 2 = 1), fuse to produce the heavy object 
(Higgs boson and diquarks in the original model). The corresponding cross-sections for qq and Higgs 
boson production read: 

<KT(s) = (-Jg-V V) I + J II {^rf 2 M%^exp(-2A JJ «?)}^T; 

V QQ / \ . — 1 9 / *' — 1 9 



««' V=l,2 
2e / 



d*H C (s) = II {^^^exp(-2A^)} 



(1) 



The variables Vi and ki respectively denote the transverse momenta of the outgoing protons and 
quarks, & are the proton fractional momentum losses 2 , and r\i are the quark rapidities, gh is the gluon- 
initiated Higgs boson production cross-section while a q ^ c is the hard qq production cross-section in the 
exclusive case. Indeed, in this exclusive process submitted to the Jz = constraint |H] 3 (a helicity 
selection rule for the production of a scalar system from the vacuum channel), the qq differential cross- 
section writes 4 , apart from normalizations included in C qq : 



exc 



i^da p(l - p) 4m£ 
99 ~ 24 dt m 2 T1 m 2 T2 ' 9 ' 1 ' 

In the model, the non-perturbative input is naturally related to the soft Pomeron trajectory taken 
from the standard Donnachie-Landshoff parametrization |1U|. namely a(t) = 1 + e + a't, with e « 0.08 
and a' w 0.25GeV- 2 . 

There are other parameters to be fixed, coming from the non-perturbative gluon propagators. Phe- 
nomenological contraints are obtained from the physical values of the total cross-sections, leading to four 
unknown parameters in formula JQl, namely the normalizations C q q, Ch and the slopes in momentum 
transfer \ q q, Ah- At this point of our study, the slopes X q q, Xh are kept as in the original papers £Q 5 . 

The problem of the normalization constants C q q, Ch requires more care, since the evaluation of 
cross-sections is the main subject of the paper. All constants of the problem (non-perturbative proton- 
gluon coupling, normalization of the non-perturbative gluon propagators, color factors, perturbative Higgs 
boson and qq production vertices) are contained in the normalizations C qq and C'h ■ The non-perturbative 
factors are poorly determined by the current data, implying large uncertainties on production cross-section 

2 The £i are given by ^i^2y/s = rrvr\ expirji + rrtT2 exp±?72, mxi = + and restricted to be smaller than 10%. 
This value for £ < 0.1 comes from the usual cut on rapidity gap to be > 2. 

3 Note that the exclusive model of formula (1) is related to the soft pomeron exchange with a low value of the intercept, 
by contrast with the model of Ref. '.I . 

4 Here, as in the rest of the paper, the hard differential cross-sections a are normalized to the usual da/dt expressions 
through 

5 They are constrained by using (an approximate parametrization of) the nucleon form factor as the Pomeron coupling 
to proton 
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"Inclusive" 

Figure 1: Double diffractive vs. standard gluon fusion production schemes, xf are the momentum fractions 
of the fusing gluons, ki are the transverse 2-momenta of the outgoing difermion in the central rapidity 
region. The "Standard" Higgs boson production is displayed for reference. "Exclusive" and "inclusive" 
double diffraction are represented in the framework of the non factorizable Pomeron model (see formulae 
(111 1171) and text for the complete kinematical notations). The hatched region represents the diffractive 
interactions at both (anti)proton vertices, while the vertical thin line is for the soft gluon exchange in the 
model. Note that xf = 1 in the "exclusive" case. 
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predictions. In particular, even if every other parameter is fixed, the non-perturbative proton-gluon 
coupling G remains undetermined and is arbitrarily fixed to G 2 /At: = 1 in the original publications P^. 
However, and it is an important aspect of our model studies, the ratio Cn/Cqq is well defined, and 
independent of everything else than the color structure of the process. This means that the ratio of Higgs 
boson to diquark production is fixed by this factor. Given the expressions of C'h and C q q pQ, one finds : 



where enter only in © the Fermi constant Qp, and the ratio of color factors to produce either a color 
singlet qq or the scalar Higgs boson in the non-factorizable Pomeron model (see for details, the second 
reference of pP). This feature will remain valid in our model of inclusive production, suggesting that the 
known, large cross-section dijet production process can be used to calibrate Higgs boson production. 

Also note that, as can be seen in (J2J, the exclusive production rate for a given quark flavour is 
proportional to its mass squared, so that light quark production is expected to be negligible, reflecting 
in another way the Jz — constraint 

2.2 Inclusive production 

The inclusive mechanism is described in the third graph of Fig. ^ The idea is to take into account that 
a Pomeron is a composite system, made itself from quarks and gluons. In our model, we thus apply the 
concept of Pomeron structure functions to compute the inclusive diffractive Higgs boson cross-section. 
The HI measurement of the diffractive structure function and the corresponding quark and gluon 
densities are used for this purpose. This implies the existence of Pomeron remnants and QCD radiation, 
as is the case for the proton. This assumption comes from QCD factorisation of hard processes. However, 
and this is also an important issue, we do not assume Regge factorisation at the proton vertices, i.e. we 
do not use the HI Pomeron flux factors in the proton or antiproton. 

Regge factorisation is known to be violated between HERA and the Tevatron. Moreover, we want 
to use the same physical idea as in the exclusive model pQ, namely that a non perturbative gluon 
exchange describes the soft interaction between the incident particles, as in Fig. ^ In practice, the 
Regge factorisation breaking appears in three ways in our model: 

i) We keep as in the original model of Ref pQ the soft Pomeron trajectory with an intercept value of 



ii) We normalize our predictions to the CDF Run I measurements, allowing for factorisation breaking 
of the Pomeron flux factors in the normalisation between the HERA and hadron colliders 6 . 

iii) The color factor derives from the non-factorizable character of the model, since it stems from 
the gluon exchange between the incident hadrons. We will see later the difference between this and the 
factorizable case. 

The formulae for the inclusive production processes considered here follow. We have, for dijet pro- 
duction 7 , considering only the dominant gluon-initiated hard processes: 

6 Indeed, recent results from a QCD fit to the diffractive structure function in HI II 21 show that the discrepancy between 
the gluonic content of the Pomeron at HERA and Tevatron 1131 appears mainly in normalisation. 
7 We call "dijets" the produced quark and gluon pairs. 




(3) 



1.08. 
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^ JJ ' V i=l,2 J i=l,2 

x{ajjGp(x 9 1 ,^G P (x 9 2 ^)}; (4) 

and for Higgs boson production: 



da ™i = Ch (^fs ^6-^) II {gp(^?,m) dxfd\ exp (-2«?A H )} ; (5) 

For the two following processes, the quark-initiated contribution can not be ignored. We have, for dilepton 
production: 



da™ 1 = C u (^r) v q ^uSW ( Vi+ki J J] {Qp(xf,M)^^fd 2 Ul d 2 fc l 6 2Q '^exp (-2«?A„)} ; (6) 
and for diphoton production: 



,1=1,2 i=l,2 



do™ 1 = C n (^-) 2 5 {2) ( Vi+ki) J] {dtithtfvicpkit^exp (-2^ 2 A 77 )} x 

^ 77 / \i=l,2 J i=l,2 

x {a gg ^ 11 Gp{x 9 1 , n)Gp{xl , fi) + a q q^ 11 Qp{x{, fj,)Q P (x?,, n)} . (7) 

In the above, xf are the Pomeron's momentum fractions carried by the gluons or quarks involved in the 
hard process, and the Gp (resp. Qp) are the Pomeron gluon (resp. quark) energy densities, i.e. the 
parton density multiplied by x\ . We use as parametrizations of the Pomeron structure functions the fits 
to the diffractive HERA data performed in |14| . The dijet cross-section 8 is now (summing over quark 
flavours /, and now including the contribution from gluon jets): 



I /Uma I 9\ f 4m Tl m T2 g ^VT\VT2 

°JJ = 2^ a ss-+q?(P ) + 108 Oga-^aaW) \ P = —Tp i P = —772 — i 

J J J J J 

„ pj d e!\ d g p f ) . J A p 9 

2 f 2 \ 9/1 1 fi I ' 99^99 — n 2 „2 



f T1 m f T2 V 2 / V 16 / Pt\Pt2 V 4 



(8) 



to be compared with 0. The above formulae are derived using |16| . and the dilepton and diphoton 
cross-sections are taken from The expressions for a gg ^ qg , <J gg ^ gg , a qq ^u, and er 5g ^ 77 , a qq ^ 11 in 

terms of the Mandelstam variables are recalled in the Appendix. 

8 The formulae JSJ are corrected for a factor 2 error coming from a known misprint in the normalization of o $ t in 1151 . 
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In the inclusive case, contrary to the exclusive case, dijet production is flavour democratic and thus 
the in JSJ extends over all flavors except for the too massive top quark, due to kinematics. Note that 
the non perturbative parameters are kept the same as in the exclusive case. Indeed, the expressions |JT]| 
can be recovered in the limit Gp — > 6(x— 1), and substituting back equation instead of © for the hard 
cross-section, which restricts to the Jz = component of the qq cross-section and reintroduces the flavor 
mass hierarchy. The normalization of Gp is not determined by the HERA data (since it is mixed with 
the flux factors) but is fixed in the model in order to lead to the same result for the energy-momentum 
sum rule as for the exclusive case. It is interesting to note that the comparison with the observed dijet 
production rates will give the correct order of magnitude for the inclusive model, while the exclusive one 
leads to too important rates 9 , as discussed in the next section. Note that the normalization cancels in 
the ratio ( 5 s - ) , with the same value as in the exclusive case, see OH 



The CDF measurement of double diffractive dijet production [5] is used as a verification of the validity 
of our approach (namely concerning the inclusive picture we consider, and the application of structure 
functions measured in electron-proton collisions in our context). Once the model validity tested, the 
measured cross-section will allow us to fix (within experimental errors) the absolute normalization of the 
cross-sections. 

In the CDF measurement, an outgoing antiproton is measured on one side of the detector, and the 
DPE nature of the events is ensured by requiring a rapidity gap on the opposite side. The selection then 
requires at least two jets satisfying a transverse energy criterion. The details of the selections can be 
found in pj. The measured cross-section is 43.6 ± 4.4(stat) ± 21.6(syst) nb, with large error bars. 

Reproducing the experimental selections on the cross-section estimates, and keeping fixed all param- 
eters as in the exclusive case we obtain a raw prediction of 11.4 nb, i.e. a factor 3.8 smaller than the 
measured mean value 10 . This prediction includes a fast simulation of detector resolution effects, using 
SHW [T§|. Considering the large uncertainties, this result is quite encouraging. Indeed, as mentioned 
above, the experimental errors are yet quite large. 

In order to verify the dynamics of the model, it is interesting to consider the dijet mass fraction, defined 
as the ratio of the mass measured in the central detector to the missing mass to the outgoing proton and 
antiproton. For exclusive events, this ratio is expected to be about 0.8 [3] given detector inefficiencies; if 
inclusive events dominate, and the model is correct, one should observe a broad distribution below one, 
essentially given by the product of the Pomeron structure functions. Fig. displays the mass fraction as 
measured by CDF, and the prediction of the present model 11 . 

Reasonable agreement is observed, suggesting that the HERA Pomeron structure functions allow for 
a correct description of inclusive DPE events. It is worthwhile to emphasize the strong influence of the 

9 One may think naively that this overestimation using parameters from soft hadronic cross-sections might be an argument 
against the non perturbative gluon model. This argument does not hold against the inclusive predictions which gives already 
the correct order of magnitude. 

10 At a theoretical level, it is interesting to note that a mere reduction of the non-perturbative proton-gluon coupling 
(arbitrarily fixed in original papers 0) from G 2 /4n = 1 — •> 1/2 swallows the normalization factor. 

11 Recall that in our predictions, we only included the gluon structure functions in the pomeron (Gp), neglecting the 
quark structure function. This assumption is justified by the weakness of the quark structure function (10 %), and the large 
error bar on the gluon structure function of about 50% 1141 . 




3 Model predictions for DPE 



3.1 DPE Dijets at the Tevatron Run I 
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(3) 


(4) 


Dijets 


7.010 5 


3.210 3 


1.210 7 


2.510 b 


Diquarks 


3.810 3 


4.4 


3.610 4 


1.710 3 


bb 


7.610 2 


8.710- 1 


9.410 3 


4.510 2 


77 


3.510- 2 


4.910" 5 


4.410" 1 


1.910~ 2 


l+l- 


4.610" 2 


6.810" 5 


1.1 


1.310" 2 



Table 1: DPE cross-section (pb) for dijets, diquarks, bb, diphotons and dileptons (fb). (1): at the 
Tevatron, \y\ < 5, m > 40 GeV; (2): at the Tevatron, \y\ < 5, m > 120 GeV; (3): at the LHC, \y\ < 5, 
m > 40 GeV; (4): at the LHC, \y\ <5,m> 120 GeV. 



behaviour of the (hard) gluon distribution of the Pomeron on the predicted mass fraction. In Fig. we 
consider a "proton-like" (i.e, soft) gluon distribution in the Pomeron, leading to unsatisfactory results. 

3.2 Model predictions for Tevatron and LHC: dijets, diphotons and dileptons 

Assuming a global normalization as determined here and all other model parameters as given above, one 
is now in situation to give predictions for Higgs boson production at the Tevatron and the LHC. In a first 
stage, we assume no evolution of this normalization with energy. We will discuss later on the influence of 
the various parameters on the predictions, analyze the (large) uncertainties which still affect the results 
and, above all, what can be done in the future to handle them, both on experimental and theoretical 
grounds. For this sake, after scaling our model to the CDF Tevatron run I predictions, we give predictions 
for various DPE processes which are relevant for Tevatron Run II and LHC measurements. 

In Figs. 0] and we give the integrated cross sections for dijets, quarks, bb, dileptons and diphotons 
at generator level for a mass higher than the mass m given on the abscissa for the Tevatron and the LHC. 
Each figure is displayed for two different mass ranges, namely between 10 and 100 GeV, and 100 and 160 
GeV. 

The differential cross-sections for dijet, dilepton and diphoton production are also given in Fig. Eland 
The dijet cross-section is dominated by the gluon contribution. The quark jet contribution amounts 
1% at small masses, and goes down to 0.1% at higher masses, both at the Tevatron and the LHC. The 
bb contribution, which is about 20% of the total quark contribution, represents only about 0.02% of 
the dijet yield. The diphoton cross-section is small (4 orders of magnitude below the bb cross-section 
both at Tevatron and LHC) but still measurable at Tevatron at low masses to test the model. These 
processes are much smaller than dijet production (due to the weak QED coupling constant, and to the 
small quark component of the Pomeron which initiates these processes), but they do have appreciable 
advantages which will appear in the following discussion. The dilepton cross-section is small (same order 
of magnitude as the diphoton) but enhanced by the presence of the Z pole (see Fig. 01 . Numerical values 
for the cross-sections are also given in Table EOl 

3.3 Model predictions for the Tevatron and the LHC: Higgs boson production 

Our results for Higgs boson production at the Tevatron and the LHC are displayed in Figs. 151 1§1 and ITU1 
and in Tables I3~3l and I3~3l In Fig. [SJ we see that the Higgs production cross-section at the LHC is higher 
by more than two orders of magnitude than at the Tevatron. In the same figure is also displayed for 
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Figure 2: Observed dijet mass fraction (CDF Run I), compared to our model prediction, using the 
Pomeron structure functions. The dashed histograms are for exclusive production with and without (c/. 
the peak at 0.8) detector simulation. 




Figure 3: Dijet mass fraction obtained using the gluon structure function from the proton. 
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Figure 4: Dijet, diquark, bb, dilepton, and diphoton cross-sections (pb) at the Tevatron. The cross- 
sections are given for a mass m above the value on the abscissa, for two different mass ranges. For 
comparison, we also display the Higgs boson cross-section for Mniggs = m- 
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Figure 5: Dijet, diquark, bb, dilepton, and diphoton cross-sections (pb) at the LHC. The cross-sections 
are given for a mass m above the value on the abscissa, for two different mass ranges. For comparison, 
we also display the Higgs boson cross-section for Mmggs = m - 
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Figure 6: Differential dijet production cross-section (pb) at the Tevatron and the LHC. The transverse 
energy of the central jets satisfies Et > 10 GeV, and their rapidity is limited to \y\ < 4. 




Figure 7: Differential diphoton and dilepton production cross-sections (fb) at the Tevatron and the LHC. 
The dilepton cross-section corresponds to a single lepton flavour. The transverse energy of the central 
particles satisfies Et > 10 GeV, and their rapidity is limited to \y\ < 4. 
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M H 


(1) 


(2) 


(3) 


(4) 


100 


3.8 


3.2 


0.3 


0.0 


110 


2.3 


1.8 


0.2 


0.0 


120 


1.3 


0.9 


0.1 


0.1 


130 


0.7 


0.4 


0.0 


0.2 


140 


0.3 


0.1 


0.0 


0.2 



Table 2: DPE Higgs production cross-section at the Tevatron (fb). (1): total cross-section, (2): bb 
channel, (3): rr channel, (4): W + W~ channel. 



M H 


(1) 


(2) 


(3) 


(4) 


100 


182.3 


152.1 


12.4 


1.5 


110 


172.6 


138.2 


11.4 


6.4 


120 


158.5 


114.3 


9.6 


18.1 


130 


147.0 


85.2 


7.2 


37.6 


140 


137.7 


54.3 


4.6 


61.6 


150 


127.5 


26.8 


2.3 


83.4 


160 


122.5 


6.2 


0.5 


109.0 


170 


115.3 


1.3 


0.1 


110.8 


180 


108.9 


0.8 


0.1 


101.4 


190 


103.8 


0.5 


0.0 


81.3 


200 


98.1 


0.3 


0.0 


72.5 



Table 3: DPE Higgs production cross-section at the LHC (fb). (1): total cross-section, (2): bb channel, 
(3): rr channel, (4): W + W~ channel. 



comparison the cross-section for standard Higgs production, which is more than two orders of magnitude 
larger. In Figs. l9l and 1 101 we give the Higgs boson production cross section at the Tevatron and the LHC 
for the different decay modes of the Higgs boson. 

The diffractive Higgs production at the LHC is mainly interesting for the lower Higgs boson masses. 
When the Higgs boson mass is heavy enough, the WW and ZZ^ decay modes become dominant and 
the visibility of the corresponding channels are already very good in standard non-diffractive events ; 
hence we do not expect the contribution from diffractive channels to be as important there. 

At lower Higgs masses, the standard non-diffractive searches are done using the 77 decay mode, which 
is loop-mediated and has very small branching fraction. In the present case, the good mass resolution 
obtained using roman pots or microstation detectors [21] allows to use the Higgs boson decays into b- 
quarks and into r-leptons, which are the two main decay modes when the Higgs boson is lighter than 
~ lAOGeV/c 2 . We will discuss this point in more detail later in this paper. 
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Figure 8: Higgs boson production cross-section at the LHC and the Tevatron. The upper plot gives 
the cross-sections at the Tevatron and the LHC for low Higgs mass to show the difference between both 
accelerators. The bottom plots show the same distribution at higher Higgs masses (the standard Higgs 
production cross-section is shown for comparison). 
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Figure 9: Higgs boson production cross-section at the Tevatron. Various decay channels are plotted as 
a function of the Higgs boson mass. 




Figure 10: Higgs boson production cross-section at the LHC. Various decay channels are plotted as a 
function of the Higgs boson mass. 
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4 Parameter dependence of DPE cross-sections 



In this section, we study the dependence of our predictions for dijet, diphoton, dilepton and Higgs boson 
cross-sections on the different parameters of the model. If we refer to the formulae given in the first 
section, the relevant parameters of our models are the momentum transfer slopes Xi (where i stands for 
H, JJ, 11, qq, 77), and e and a' respectively the intercept and slope of the Pomeron Regge trajectory. 
We will discuss the dependence on e in the next section (note that this value is however imposed in our 
model, but changes in other Pomeron models, as discussed later on). 

The dependence on Xh, A; (where i stands for 11, J J, bb, qq, and 77) is displayed in Tables 0] (for 
Tevatron) and 0] (for LHC). In our model, we assumed Xh = 2 and A, = 3 (as in the original exclusive 
model 1 ). We now vary these values by one unity, keeping values (assuming or not the equality of 
the A's) around those physically connected to the nucleon form factors. We also give numbers with the 
pomeron slope a' taken to be 0.1. We note as expected that the cross-sections for dilepton, diphoton, 
dijet and bb production vary by the same factor when we vary the parameters. 

The results of Tables 0] and 0] are shown for each value of the A's (except the reference one on the 
first row) on two lines: on the first line, we give the cross-section values as they come directly from the 
generator simulation, and on the second lines we rescale the values to the reference dijet cross-section 12 . 
This method allows us to determine the error on the Higgs production cross-section due to the assumptions 
made on the respective values of Xh and Aj. For the Tevatron, the cross-section we obtain for a Higgs 
boson mass of 120 GeV is 1.2 fb, and varies between 0.6 and 1.9, which gives an incertitude of about 50% 
on the Higgs boson cross-section prediction. For the LHC, the default value is 160 fb and varies between 
89 and 220, which gives also an error of about 50%. 

The other assumption in our study is to take the gluon density in the pomeron at Q 2 = 75 GeV 2 , 
since this is the upper value of Q 2 for the QCD fit. If we take a value at a higher Q 2 by using a DGLAP 
evolution we obtain that the gluon density can vary by a factor 2 at Q 2 = 4000 GeV 2 . This has a clear 
effect on our predictions which can thus vary by a factor 2. 

Another assumption of our model which needs to be tested against data is that the fraction of events 
with tagged protons is the same at the LHC and the Tevatron, which in other words means that the 
gap survival probability does not depend much on the center of mass energy. To verify this point will be 
possible only when the LHC energy range will be available to experimentation. 

To summarize, we estimate that, on the basis of our model calculations, the errors due to the parameter 
dependence on our predictions for Higgs cross section (once normalized to dijets) are of the order of a 
factor 4. 

5 Testing and comparing models at the Tevatron Run II 

Let us discuss and propose ways to verify our model and determine its parameters more precisely using 
the Tevatron run II data for dijets, diphotons and dileptons. This will allow to make improved predictions 
for e.g. Higgs boson production later on at the LHC. In the same time, this will give the possibility to 
differentiate, and thus compare the existing models which give predictions for diffractive Higgs production. 

12 We remind that we scale our cross-section predictions to the CDF RunI measurement. 
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Param. 


&dijets 


a bl 






G Higgs 


Xh = 2, \i = 3 


7.2 10 B 


2.5 10 4 


1.6 


1.6 10 2 


1.2 


X H = 3, A, = 4 


4.7 10 6 
7.2 10 6 


1.6 10 4 
2.4 10 4 


1.0 
1.5 


1.2 10 2 
1.8 10 2 


0.7 

1.1 


Ah = 2, Aj = 2 


1.3 10 Y 
7.2 10 6 


4.5 10 4 
2.5 10 4 


2.8 
1.5 


3.3 10 2 
1.8 10 2 


1.2 
0.7 


A ff = 4, Aj = 4 


4.7 10 a 
7.2 10 6 


1.6 10 4 
2.4 10 4 


1.0 
1.5 


1.2 10 2 
1.8 10 2 


0.42 
0.64 


Ah = 1, Aj = 2 


1.3 10 Y 
7.2 10 6 


4.5 10 4 
2.5 10 4 


2.8 
1.5 


3.3 10 2 
1.8 10 2 


3.5 
1.9 


a' = 0.1 


9.8 10 B 
7.2 10 6 


3.4 10 4 

2.5 10 4 


2.2 
1.6 


2.5 10 2 
1.8 10 2 


1.8 
1.3 



Tabic 4: Numerical values of the cross-sections (fb) for different values of the parameters, at the Tevatron 
(the subscript i stands for 11, bb, J J, 77). The Higgs cross-section is for a Higgs mass of 120 GeV. All 
other cross-sections are given for a dijet, dilcpton, or diphoton mass greater than 10 GeV. The first lines 
read the direct output of the program, whereas the second lines are rescaled to the same value of the 
dijet cross-section, namely the CDF measurement. 



Param. 


O 'dijet s 


a bb 




era 


&Higgs 


Xh = 2, Xi = 3 


4.2 10 Y 


1.1 10 b 


8.8 10 3 


1.6 10 3 


1.6 10 2 


X H = 3, Ai = 4 


2.7 10 7 
4.2 10 7 


7.5 10 4 
1.2 10 5 


5.5 10 3 

8.6 10 3 


1.0 10 3 
1.6 10 3 


9.0 10 1 
1.4 10 2 


Xh = 2, Xi = 2 


7.1 10 Y 

4.2 10 7 


1.9 10 5 
1.1 10 5 


1.6 10 4 
9.5 10 3 


2.8 10 3 
1.7 10 3 


1.6 10 2 
9.5 10 1 


X H = 4, A 4 = 4 


2.7 10 7 
4.2 10 7 


7.5 10 4 
1.2 10 5 


5.5 10 3 

8.6 10 3 


1.0 10 3 
1.6 10 3 


5.7 10 1 
8.9 10 1 


X H = 1, ^ = 2 


7.1 10 Y 

4.2 10 7 


1.9 10 5 
1.1 10 5 


1.6 10 4 
9.5 10 3 


2.8 10 3 
1.7 10 3 


3.8 10 2 
2.2 10 2 


a' = 0.1 


6.1 10 Y 

4.2 10 7 


1.7 10 b 
1.2 10 5 


1.5 10 4 
1.0 10 4 


2.3 10 3 
1.6 10 3 


2.5 10 2 
1.7 10 2 



Table 5: Numerical values of the cross-sections (fb) for different values of the parameters at the LHC (i 
means 11, bb, J J, 77). The Higgs cross-section is for a Higgs mass of 120 GeV. All other cross-sections are 
given for a dijet, dilcpton, or diphoton mass greater than 10 GeV. The first lines read the direct output 
of the program, whereas the second lines are rescaled to the same value of the dijet cross-section, namely 
the CDF measurement. 
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5.1 Brief discussion of other models 



Models attempting to describe central diffractive production from pp or pp interactions base their predic- 
tions on either an explicit color singlet exchange of two gluons (where one may be soft, as in our model), 
or in terms of hadronic interactions, in which diffractive features (rapidity gaps, leading protons) appear 
in relation with soft initial and/or final state interactions. 

In this section we attempt to summarize results obtained in different pictures for the same observables, 
and assess how observation could allow to distinguish them. We concentrate on recent models 13 which 
have provided explicit numbers for various processes of interest (namely dijet, diphoton and Higgs boson 
production) so that a comparison of predictions is made possible. We therefore compare our results with 
the exclusive model predictions of Ref. |22|-a. with the double Pomeron inclusive model of Refs. |22|-b 
and with the soft color interaction models of Ref. |22J-c. These three approaches will be shortly described 
below. 

Two-gluon models for exclusive production have been first proposed, (c/., see 2.), but led in general 
to too large cross-sections, in conflict with the upper bound from CDF dijets. Despite their initiatory 
role in the problem, they will not be discussed in detail here. We cannot either include the original 
Bialas-Landshoff exclusive model .1 in this discussion, since the model has not been extended to gluon 
pair production. A confrontation of the exclusive Bialas-Landshoff dijet production with the CDF upper 
bound is thus not possible at this point. 

Exclusive production is evaluated in |22j-a as a perturbative process involving two-gluon exchange 
tamed by Sudakov suppression factors supplemented by rapidity-gap suppression. Two gluons are ex- 
tracted from the beam particles according to the usual proton structure functions (the process can thus 
be interpreted as "proton-induced") and one of them couples to the hard central system (e.g. a high 
mass dijet, Higgs boson etc.). The infrared divergence of this process (da cx dQ^/Qf, Qt being the 
gluon transverse momentum) is regularized by the requirement that no radiation occurs along the gluon 
which would fill the rapidity gap. This leads to Sudakov-like form factors effectively damping the \ow-Q t 
region, so that the whole process can be considered as entirely perturbative and thus calculable. Finally 
a conventional rapidity gap survival factor corrects for the soft rescattering corrections. Clues for future 
studies on distinguishing the exclusive from the inclusive DPE processes are given in Section 15.21 

The model for inclusive DPE presented in [22|-b, and implemented in PDMWIG 23 is a direct transpo- 
sition of the Pomeron flux and parton densities from the ep to the pp context, and is therefore referred to 
as a factorizable model 14 . In particular, it uses the full Regge parametrization of Pomeron flux factors as 
determined on the HERA data quoted before. The factorizable model is formally similar to ours, since 
both models satisfy QCD factorization by the use of gluon structure functions in the Pomeron to describe 
the hard process in the central rapidity region. However, the choice of parameters and color structure 
reflects a different interpretation of the nature of the Pomeron, as outlined in section 15.31 The closely 
related formulation of both models will allow to discuss their differences in terms of physical parameters. 

Two models have been designed to describe diffractive physics as soft color rescattering over a hard 
subprocess [22J-C, namely SCI (Soft Color Interaction) and GAL (General Area Law). They are imple- 
mented as a transition between the hard interaction and hadronization, and therefore fit naturally in 
Monte-Carlo programs such as PYTHIA [23]. In these models, color exchanges occur in the final state, 

13 In fact, the observable tests that we propose could be applied to other models as well, and thus have a wide applicability. 

14 As in our model, the model is Pomeron-induced. However, our model is based on a soft Pomeron interaction between 
the initial hadrons which do not verify Regge factorization at the incident vertices. Here the model is based on a Regge- 
factorized mechanism with a Pomeron at each vertex, the necessary factorization breaking coming from the gap suppression 
factor. 
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potentially stopping color flow between the remnants of the incoming hadrons and the central system, 
leading to diffractive event topologies. These final state interaction models have a different formulation 
from the previous ones 15 . We will restrict ourselves to recall their predictions for various DPE final states, 
insisting on differences with the other models. 

5.2 Exclusive "proton-induced" model 

The question of the observation of exclusive DPE events is open. As we have said, the current data 
allow only to place an upper bound on the process, and this bound is still compatible with (but not far 
lower from) the predictions of 22 -a. The forthcoming Tevatron Run-II data will give new insight on 
this question, by looking for diffractively produced central dijets, or even better diphotons and dileptons 
since the events will be cleaner. 

A specific question concerns the separability of the exclusive process from the inclusive background. 
The exclusive processes will show a high value of the dijet, dilepton or diphoton mass fraction. In our 
model, these events correspond to the case when the pomeron remnants show very little energy. It is thus 
worthwhile to verify whether one is able to distinguish between these two configurations. If we assume 
an experimental resolution of 15% on the measured dijet mass, one should require Mf rac > 0.85 for both 
the exclusive or the inclusive events. With this cut, we obtain a cross-section from our inclusive model of 
about 2 fb at LHC, which is large enough to be competitive with the exclusive numbers given in |22j-a. 
The Higgs boson cross-section as a function of the mass fraction is given in Fig. El hi fact, this shows 
that it will be difficult to distinguish between pure exclusive events and the tail of the inclusive events 
which we can call "quasi-exclusive" events 16 . It is thus important to be able to tag both kinds of events 

ED 

Another example of discriminating observable is given by the ratio of the diphoton and dilepton cross- 
sections, as a function of the mass-fraction value. In the inclusive models, this ratio is determined by 
the quark and gluon distributions inside the Pomeron, and the presence of the Z pole in the dilepton 
cross-section. It is illustrated in Fig. The presence of an exclusive contribution to the cross-section 
would result in a sharp enhacement at Mf rac ~ 1, since diphotons are allowed in the exclusive case, 
whereas dileptons are absent. 

In any case, measurements of the cross-section of exclusive or quasi-exclusive events at Tevatron, in 
the case of dilepton, diphoton or dijet production will provide a direct test of the models. 

5.3 Comparison of the "Pomeron-induced" models at the Tevatron Run II 

It has already been noticed that our model presented here, and the factorizable Pomeron induced model 
of [22]-b in practice differ in the relative normalizations of Higgs boson vs. dijets and in the value of the 
parameters. In particular, we take e = 0.08, see formula 0-0 while the factorizable model assumes 
a value of the pomeron intercept e = 0.2 coming from the HERA measurements. Also, the slope in 
momentum transfer originates from the nuclcon form factor in our model (eventually supplemented by 
an additional term related to the non-perturbative gluon propagator) whereas the factorizable Pomeron 
induced model takes the value from HERA measurements. It is thus interesting to vary this parameter. 

15 Note a formal similarity with the Bialas-LandshofT mechanism with the superposition of hard and soft color exchanges 
in the same probability amplitude. 

16 Let us mention another source of backgrounds to the exclusive process, coming from the experimental difficulty of 
telling that a tagged forward proton is not resulting from an N* decay and accompanied by other undetected particles. 
Provided the decay proton is still in the required momentum range (£ < 0.1), this will result in an underestimation of the 
missing mass, i.e. an overestimation of the mass fraction, hence polluting the selection. 
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Figure 11: Remaining diffractive Higgs boson cross-section at the LHC, as a function of the cut on the 
mass fraction. 
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Figure 12: Diphoton to dilepton cross-section ratio, as a function of the mass fraction. 
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Note an important difference between model predictions due to the distinctive color reconnection 
pattern of the two models. In the factorized model the Higgs boson to dijet normalization ratio is: 
( £b- ] = v ?. Gf , with a factor 8 smaller than O) for the non- factorized model. The reason is that the 

\ CjJ Jfact 2477 

factorizable model authorizes dijets coming from two quarks or gluons bearing different color charges, 
while they are restricted to be in an overall singlet state in the non-factorizable model due to the com- 
mon gluon exchanged between the incident nucleons (in other terms the whole of Pomeron remnants is 
essentially color singlet). Hence the prediction of the Higgs boson vs. dijet cross-section ratio keeps track 
of this difference. Here again, insight can be gained from the comparison of diphoton (or dilepton) and 
dijet cross-sections: since the first ones have the same color structure as Higgs boson production, they 
can be used to infer which, of the factorizable and non-factorizable pictures, is correct. 

Fig. shows the variation of the Higgs and dijet cross-sections as a function of e at the generator 
level, and their ratio. As can be guessed from formulae Q-Q, an increase of e will enhance the cross- 
section essentially at low diffractive mass, therefore the effect on the dijet cross-section is much larger 
than on the Higgs boson cross-section. When e is increased, the dijet cross-section is increased, and thus 
the multiplicative factor needed to adjust the prediction to the CDF Run I measurement is smaller. This 
explains why the predicted Higgs cross-section is expected to decrease when e increases. 

Fig. 1141 shows the e-dependence of the differential dijet cross-section dajj/dMjj, when compared to 
(i.e., divided by) a reference distribution, taken at e — 0.08. A strong dependence is found, meaning that 
the study of this distribution, independently from its normalisation, should allow to constrain the value of 
e, and to obtain a more precise prediction for the Higgs boson cross-section at the LHC 17 Similar studies 
can be done with dilepton and diphoton events, which will benefit from much central mass resolution. 

In Fig. El we give the cross-sections for dijets, bb, 77, and Higgs boson, and in Fig. El the ratio of 
the Higgs cross-section with respect to 77 and bb as a function of e for the LHC. Namely, it will allow to 
get a prediction on these cross-sections, once e been determined, e.g. at the Tevatron. 

These last remarks highlight once more the interest of DPE diphoton and dilepton production. In- 
deed, the experimental resolution on photon, electron and muon momenta is much better than the jet 
energy resolution, since in this case the events benefit from tracking information and/or electromagnetic 
calorimetry with a much better defined energy scale. Moreover, especially in the dilepton case, the hard 
sub- process is essentially quark-initiated (see formula iJSJ), and the knowledge of Qp is much more pre- 
cise than that of Gp, thus eliminating another source of uncertainty. The measurement of e will also 
be possible using the diphoton and dilepton measurement and the dependence of the cross-sections as a 
function the diphoton or the dilepton mass. 



5.4 Predictions: summary and comparison 

In Table lBT4l we give, where available, the cross-sections obtained by the different models. The factorizable 
numbers at the LHC (second column) should contain a factor accounting for the unknown rapidity gap 
survival probability (denoted SP); this factor is around 0.1 at the Tevatron. The numbers from the 
proton-induced model refer to exclusive production, and are as such not refering to strictly the same 
process as the pomeron induced ones. The various models have fundamentally different philosophies and 
assumptions, so that caution should be taken not to conclude from any numerical agreement that the 
models themselves agree. 

17 The sensitivity of dcrjj/dMjj to e will however crucially depend on the experimental determination of the jet energy 
scale and on the remaining uncertainties on the Pomeron structure function Gp. 
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Figure 14: Differential dijet cross-sections for different values of e, compared to a reference taken at 
e = 0.08. All cross-sections are normalised to the same value. The values e = 0.08, 0.2 and 0.5 correspond, 
respectively, to the soft, hard (measurement at HERA), and hard BFKL pomeron intercepts. 
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Figure 15: Raw cross-sections vs. e 
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Figure 16: Higgs compared to dijets and diphotons as a function of e 
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process and cuts 


(1) 


(2) 


(3) 


(4) 


Higgs, 115 GeV, Tev 


1.7 


0.029-0.092 


0.03 


0.00012 


Higgs, 115 GeV, LHC 


169. 


379.-486. xSP 


1.4 


0.19 


Higgs, 160 GeV, LHC 


123. 


145. xSP 


0.55 




77, Tev, E T > 12GeV, r) < 2 


71. 


128. (27.) 




- 20. 


77, Tev, E T > 12GeV, r] < 1 


9. 


8. (2.) 






77, LHC, E T > 50GeV, r) < 2 


1.5 






0.1 


77, LHC, E T > UOGeV, rj<5 


19. 




0.12 





Table 6: DPE Higgs and diphoton production cross-section at the LHC and the Tevatron (fb) for different 
models. (1): non factorisable pomeron based model (present work), (2): factorisable pomeron based 
model (the results contain both the pomeron and the reggeon components. For the 77 cross-section, we 
give in parenthesis the contribution of the pomeron only). (3): proton based model (exclusive case; the 
results are given for a Higgs mass of 120 GeV), (4): soft color interaction models. 

6 Advantage of diffractive vs. "standard" Higgs boson produc- 
tion 

As is now well-known [2 El El (HJ-a, exclusive DPE events are very attractive from the experimental 
point of view, since the measurement of the outgoing protons momenta, using dedicated forward detectors 
(like roman pots or microchambers [2]) allows to determine the mass of the centrally produced system 
with excellent precision, so that the Higgs boson signal can be extracted from the DPE continuum back- 
ground. Moreover, the background itself is strongly suppressed due to angular momentum conservation 
constraints (the u Jz = 0" rule [5], see section l2~ll . Exclusive DPE events are however not an experimen- 
tally established process according to the data currently at our disposal. As we have emphazised, the 
"quasi-exclusive processes" which are quite close in topology from the exclusive processes are established 
processes and may keep an essential part of the nice kinematical features of the exclusive production. 

In general inclusive DPE events, the angular momentum constraint above does not apply and the 
background is not suppressed, so that the a priori signal-to-background ratio (S/B) is similar to that of 
"standard" non-diffractive events. The missing mass method will not work as such, since the equality 
between the missing mass to the outgoing protons and the mass of the central system does not hold. 

However, inclusive DPE events can still be kinematically constrained in a much stronger way than 
non-diffractive events. If, in addition to the forward proton detectors, the experiments can dispose of 
very forward calorimetry (e.g. up to 77 ~ 8, as was proposed in [53]), one may use the measurement of 
the Pomeron remnants to fully constrain the events kinematics. Combined with the higher cross-sections 
predicted in this last case, it shows that such forward calorimetry would be adapted and necessary. The 
first resolutions using these detectors are given in 

The idea is to use the forward calorimeters as vetos. The central mass is reconstructed with the 
proton momenta only, and events are considered only if the measured remnant energy is small enough. In 
a sense, this is a way to select "quasi-exclusive" events. Fig. 3 of Ref. [S] shows the achievable resolution 
as a function of the remnant veto. 

Let us finish this discussion with an important point concerning the interest of the rr decay mode. 
According to section the inclusive DPE tt cross-section is very small compared to the DPE dijet 
cross-section, and this is understood from the smallness of the QED coupling constant and the quark 



24 



component of the Pomeron. Compared to this, the Higgs boson branching fraction into r pairs is still 
~ 10% in the mass range we consider. It follows from the cross-sections given earlier (see Table 13.31 
and Fig. I1(J|) that the cross-section ratio for r-pair events is 0(10) at for example m TT ~ 130 GeV for 
inclusive DPE events, whereas this mode is invisible through non-diffractive events. 

Of course, the missing energy in r decays, the presence of the Z peak in the r-pair cross-section, and 
other background processes not considered here spoil these numbers to some extent, but we do consider 
that this channel needs investigation. Also note that a pseudoscalar Higgs boson A, which appears in 
models with an extended Higgs sector (such as the supersymmetric models), does not couple to gauge 
bosons and hence preserves a significant branching fraction into r-pairs up to very high masses, where 
the effect of the Z peak is negligible; the same remark is true for the MSSM h boson in some regions of 
the parameter space. 

7 Summary 

The present work is meant to serve as a reference to the non factorizable, Pomeron induced DPE model. 

• The original exclusive DPE formulation are recalled, and the inclusive expressions are introduced, 
see formulae I|4I7|I . 

• A systematic study is performed, that sheds light on the role of the various model parameters. It 
is emphasized that the forthcoming Tevatron data will allow to constrain their values, and answer 
the pending question of the existence of exclusive DPE production. 

• Precise predictions for discovery physics at the LHC will then be possible; we focused on the Higgs 



• It is recalled that, assuming favourable exprimental installations, the bb decay mode of the Higgs 
boson may well be exploited in double diffractive production. 

• On the theoretical side, the H— > tt channel is most promising, although some experimental diffi- 
culties have to be overcome. 
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Appendix 

The cross-sections for the hard processes that are used in our study are summarized below, in terms of 
the usual Mandelstam variables. For dijet production, we take: 



boson case. 




(10) 



(9) 
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For dilepton production, we have: 



daqq-.ii 7ra 2 4 u 2 + t 2 

Here Cew includes electromagnetic and weak factors accounting for photon and Z exchange between the 
initial quark pair and the final leptons. Its expression is, for a quark of flavour i: 



2 _ s( S ~M 2 )(L e + R e )(L q + R q ) s 2 (L 2 e + R 2 e )(L 2 +R 2 q ) 

6i 8x w (l - x w )[(s - M 2 Z ) 2 + M*T%] + 644,(1 - x w ) 2 [{s - M 2 ) 2 + M%Y\\ 




Above, xw is the weak mixing angle, Mz and Tz are the Z boson mass and width, and L and R are the 
left and righthanded chiral couplings. Finally, the diphoton production cross-section reads: 

(14) 

where C; oop is a complicated expression resulting from the loop integral occurring in the gg — > 77 process. 
Its detailed expression can be found in |18j . 
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